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FOREWORD 


This final report documents the work performed at the Dynamic Science 
Division of Marshall Industries under NASA Contract Number NAS7-442, 
relating to the description of the flame structure of thermally unstable fuels. 
At Dynamic Science, this work was assigned the project number SN-81 B & C. 
Two additional reports will be published under this NASA contract which will 
describe the work performed under Dynamic Science project number SN-81 A 
and SN-81 D. 

All phases of this contract were monitored by Dr. Charles Feller and 
Dr. Richard Prlem of NASA-Lewis Research Center. At Dynamic Science, 

Dr. B. P. Breen was Program Manager of the total contract (NAS7-442). 

Mr. R. J. Hoffman acted as Technical Manager of the work reported in this 
report (SN-81 B & C). The analytical model (Appendix A) was developed by 
Mr. S. Hersh, Mr. B. R. Lawver designed and constructed the experimental 
equipment and the experiments were conducted by Mr. Hersh and Mr. Lawver. 
Dr. E. W. Nadlg was responsible for the physical properties, Appendix B. 

Thanks are expressed to Mr. John White, Mr. Wm. Irwin, and 
Mr. Frank Cummings for their assistance in constructing and operating the 
experimental equipment. 


ABSTRACT 


The results of an investigation of the hydrazine/nitrogen 
tetroxide droplet flame are reported. The results of the experi- 
mental work are the temperature and stable species concentration 
profiles. It is concluded that a hydrazine decomposition flame 
exists at the droplet surface causing a large temperature gradient 
which results in much higher burning rate than found for thermally 
stable fuels . 
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SUMMARY 


The flame structure of a 2500 micron hydrazine drop burning in 
nitrogen tetroxide was determined experimentally by measuring the tempera- 
ture and specific concentrations at a plane 90 from the stagnation point. 

The hydrazine drop diameter was maintained constant by suspending it on 
the end of a water cooled hypodermic needle attached to a constant flow 
rate syringe. ^ 20 ^ was fed into the burner, establishing a thermal convec- 
tive flow field around the drop. The flame was probed with a shielded 
0.002 inch platinum-platinum/1 0% rhodium thermocouple as well as a 
concentration sampling probe. Both an "outer" and an "inner flame" were 
visually observable. The locations of these "flames" do not coincide with 
the high temperature points in the drop combustion zone, and therefore are 
not true flames but regions of strong visible radiation. 

The experiments indicate that hydrazine decomposes very close to 
the liquid drop surface producing a steep temperature gradient which controls 
the burning rate. Although not proven conclusively, the zone between the 
drop and the inner flame appears to consist primarily of hydrazine decom- 
position products. 
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INTRODUCTION 


An understanding of the processes by which a liquid fuel drop Is 
vaporized and burned In the presence of an oxidizing media Is essential 
to predicting the combustion process within a liquid rocket engine. In the 
course of an Investigation of the rate of burning of a hydrazine drop In an 
atmosphere of N O. /NO-, It was found that conventional methods of treating 
droplet burning aid not give results In agreement with the experiment. 

In an attempt to account for the difference In burning rate of an 
exothermic fuel drop (N2H^) and an endothermic drop (Heptane) , experimental 
hydrazine decomposition rate data were applied to the prediction of combus- 
tion within the low velocity regions of a rocket engine (Ref. 1). While this 
approach showed promise of predicting correct burning rate and engine per- 
formance for the hydrazine family of fuels. It appeared appropriate, in view 
of the absence of similar studies , to investigate the two flame nature of 
hydra zine/nltrogen tetroxlde burning In detail. 

The objective of this work was thus to provide a description of the 
oxidiation of thermally unstable rocket fuels so that model of burning rate 
within a Combustion chamber could be developed. 
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SYMBOLS 


D.. 

ij 


H. 


H 


fo 


K 


K 


m 


R 


W. 


W, 


Oi 


LJ 


a 




i3, 


x,m 


€. 


th 

= the binary diffusion coefficients of i specie into 
each of the other species (cm®/ sec) 

= the enthalpy of the 1^^ species at temperature T. 

(cal/mole) 

= the enthalpy of the fuel at its initial conditions (cal/mole) 

= thermal conductivity (cal/cm sec ®C) 

th 

= equilibrium constant for the m reaction 

= index for atomic species (e.g. JI>=1 indicates the 1st 
atomic specie) 

= partial pressure of the i^^ specie (P^/P = 

= gas constant 

= total number of atomic species 
= the fuel burning rate in (moles/second) 

= rate of transport of the i^^ chemical species in (moles/second) 

= the mole fraction of the 1^^ specie (moles l/moles 

th 

= the number of i atoms in the i specie 
= the number of £< atoms in a fuel molecule. 


stoichiometric coefficients for the reactants of the m 
reaction 

stoichiometric coefficients for the products of the m 
reaction 


th 


th 


V 


= total number of chemical species 


3 



DISCUSSION 


Although semi-analytical models for the determination of the burning 
rate of thermally stable fuel drops have been developed (Ref. 2 ) , an analytical 
determination of the burning rate of a thermally unstable drop, such as hydra- 
zine has not been achieved. 

Liquid hydrazine burning in NTO exhibits two visibly discrete "flames" 
(Figure 1 ). Initial studies ascribed the "inner flame" to the energy released 
by hydrazine decomposition, the "outer flame" being the region of oxidation 
of the decomposition products. 

As part of a program to gain understanding of the processes underlying 
the hydrazine/NTO flame structure an analytical model which could be used 
to determine the N„H. burning rate under several simplifying assumptions was 
formulated (Appendix A) and experimental studies of the flame properties were 
conducted . 

The transport properties required by the analysis, consisting of both 
measured, values and approximate formulas where no experimental data was 
available, are compiled in Appendix B. This information is of a general nature 
and is therefore of use for any analysis for which information of this type is 
necessary. 

Two series of experimental tests were conducted. Both series utilize 
the same basic equipment , the difference being in the instrumentation used 
to measure the flame parameters. Temperature profiles from the ambient 
atmosphere to the drop surface were obtained using a thermocouple probe while 
Individual specie concentrations in the hydrazine-nitrogen tetroxide flame were 
determined from a mass spectrometrlc analysis of samples withdrawn from the 
flame with a sampling probe. 

Hydrazine was chosenforuse in these tests because it is the basic 
ingredient in the most widely used thermally unstable fuel and the result 
would have the greatest practical application. Burning droplet tests were 
run with propellant grade hydrazine burning the oxidizers, nitrogen tetroxide, 
oxygen, air, and oxygen/nitrogen mixtures. Temperature profiles were 
measured with all of these oxidizers, however, concentration probing was 
done only for the nitrogen tetroxide oxidizer. 
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APPARATTJS AND PROCEDURE 


ObjectlvGS » — The purpose of these experiments was to measure the 
flame structure of a hydrazine droplet burning in nearly stagnant nitrogen 
tetroxide vapor. The flame structure was defined by measuring the tempera- 
ture profile and chemical species concentration in the free convective flame 
surrounding a constant diameter hydrazine droplet. Tests were also conducted 
with hydrazine burning in oxygen/nitrogen mixtures to measure the effect of 
the oxidizer concentration on the flame structure. 

These tests were conducted with propellant grade hydrazine and nitrogen 
tetroxide. The oxygen and nitrogen gases were supplied from standard gas 
cylinders . 

The Burner . - The tests were conducted with the burner apparatus* 
shown in Figure 2. Basically the burner consists of a test section which per- 
mits observation of the combustion and the related propellant feed systems. 

The burner is designed such that either liquid fuel droplets or gaseous fuels 
can be burned in gaseous oxidizers. A flow schematic of the burner propellant 
feed system is shown in Figure 3. 

The droplet flame is produced by suspending the fuel droplet from a 
water cooled hypodermic needle in the burner test section. The fuel is fed 
from a syringe into the suspension needle with a variable speed syringe drive. 
The droplet diameter can be maintained constant by proper adjustment of the 
syringe drive. 

The fuel droplet suspension needle is shown in Figure 4. It is a 
0.010" O.D. (0.005" I.D.) tube surrounded by two larger concentric tubes 
through which water is circulated. The water provides cooling of the exposed 
proportion of the small needle, to prevent boiling and thermal decomposition 
of the fuel within the needle. The oxidizing gases are fed into the burner 
through an injector tube located in the burner base plate. The oxidizer flow 
rates are controlled with calibrated sonic orifices , by regulating the orifice 
upstream pressures. Oxygen is fed to the burner orifice through a pressure 
regulator mounted on the supply bottle. The N„0 vapor is generated by 
electrically heating the storage tanks and supply lines. The N20^ vapor 
temperature and flowrate are controlled by regulating the heater current. 

The ^20^ flowrate is also measured with a calibrated sonic orifice. 


♦This burner was designed and built under Contract AF 04(611)-! 161 6. 
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The burner test section has opposing pyrex windows for observing and 
photographing the flames. The combustion products are removed from the 
burner tube with an air driven ejector mounted to the burner top plate. 

Temperature Probe. - The temperature probe is illustrated in Figure 5. 
This probe consists of 0.002" diameter platlnum/platinum-10% rhodium thermo- 
couple wires mounted in a small ceramic insulator encased in a stainless steel 
sheath. The steel sheath protects the exposed wires from the hot combustion 
gases when the probe is fully submerged into the flame, as is illustrated in 
Figure 7. 

A 0.010" diameter platlnum/platinum-10% rhodium thermocouple without 
the steel sheath was used for the initial probing of the N_H flame. These 
tests indicated that heat conducted along the thermocouple leads from the ex- 
posed portion of the wires caused erroneous temperature measurements, there- 
fore, it was necessary to provide a steel sheath. In addition the response of 
the larger thermocouple was not satisfactory. To improve the response time an 
attempt was made to construct thermocouples with 0.001" diameter wires, how- 
ever this proved to be too time consuming to be warranted. The response of the 
0.002" diameter thermocouple is on the order of 40 milliseconds. 

The thermocouples were calibrated by immersing the probe in compounds 
with known melting point temperatures. The calibration data points are shown 
in Figure 6 along with a plot of the standard calibration for thermocouples 
obtained from Reference 3. The good agreement between the observed and 
standard calibrations encouraged the use of the latter in the higher tempera- 
ture range where calibration measurements were not made. 

The thermocouple probe is traversed through the droplet flame with a 
variable speed electric motor drive mechanism as shown schematically in 
Figure 7. The thermocouple position is indicated on an oscilloscope by using 
the signal from a 5000 O linear potentiometer to drive the oscilloscope hori- 
zontal sweep. The thermocouple output voltage is applied to the oscilloscope 
vertical amplifier input, providing direct temperature versus distance profile 
throughout the flame region. The temperature profile is recorded by photo- 
graphing the oscilloscope screen with a bezel mounted Poloroid camera. 

Concentration Probe. - The concentration sampling technique described 
in Reference 4, was used to measure the chemical species concentrations 
in the flames. The technique consists of klnetically freezing the reacting 

gases by drawing them through a very small supersonic nozzle into a vacuum 
environment. The rapid expansion quenches the reactions, thus preserving 
the stable chemical species for subsequent analysis. 
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The concentration sampling probe Is a 6mm O.D. quartz tube with 
1mm walls, drawn to a fine point at the sampling end to minimize the flame 
disturbance. The nozzle throat diameter is 0.002 inches in diameter which 
gives an expansion ratio of 6200/1. Tlie inlet area ratio is 300/1. Based on 
a flame gas velocity of 6-10 ft/sec, the nozzle residence time is approxi- 
mately 40 n sec. 


The sampling system is shown in Figure 8. The gas samples were 
collected in bottles and then analyzed on a Hitachi Perkin-Elmer mass 
spectrograph. 

The mole fractions were calculated from the relationship 



where 

x^ = mole fraction of species i 
P. = partial pressure of species i 
P^ = total pressure of sample 

The partial pressures are determined from the magnitude of the mass peaks 
measured with the mass spectrograph. The partial pressures are found from 
the relationship 


n 



where h = measured peak height of mass number (m) 

m 

= sensitivity of species (i) to mass number (m) 

There is one equation for each mass number and the problem then is the solu- 
tion of a system of simultaneous equations. The number of equations is equal 
to the number of species involved. To simplify the data reduction, these 
equations were programmed on a CDC 3600 computer. The input consists of 
the measured mass peak heights of the species of interest and the sensitivity 
factors. The output is the species mole fractions. 

Analysis of the mass peaks for the hydrazine/NTO flame indicated 
that the mass numbers of interest are. 
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Mass Number 


2 

17 

18 
28 

29 

30 


32 

°2 

44 

N^O 

46 

NO^ 


Tlie mass number 29 was used for the hydrazine to distinguish it from oxygen, 
mass number 32. The hydrazine exhibits a 29 peak that is 43% of its 32 peak, 
also it exhibits a strong peak at mass number 31, The above mass numbers 
and sp.ecles were used in determining the hydrazlne/NTO flame structures. 


Species 

«2 

NH„ 


« 2 ° 

"•a 

'' 2«4 

NO 


Experimental Results 

Temperature Profiles. - A typical temperature profile for a hydrazine- 
nitrogen tetroxide flame is shown in Figure 9. The temperature begins to rise 
before the "outer flame" is reached and rises to a peak between the "inner" 
and "outer flame. " This peak temperature corresponds to a blue region which 
is visually observable but does not appear in the flame photographs . Attempts 
to use film with a greater spectral sensitivity in the blue region than the 
Polaroid color film normally used also failed to reproduce this region. The 
temperature decreases smoothly from its maximum, through the "inner flame, " 
to a relatively high value close to the drop surface. It should be noted that 
the locations of the visible "flames" do not coincide with the high tempera- 
ture points. Thus what have been denoted as the "inner" and "outer flames" 
are not true flames but regions of strong visible radiation. 

The shape of this temperature profile, in particular the high temperature 
close to the drop surface would indicate that hydrazine decomposition had 
occurred at the droplet surface. To further examine this behavior, hydrazine 
droplets were burned in oxygen -nitrogen atmospheres at various 02~^2 
tratlons. A hydrazine-air profile is shown in Figure 10. As can be seen, both 
the maximum and drop surface temperatures, are less than that obtained for 
the hydrazine-NTO system, (Figure 9). 


8 


The 02~^2 were varied over a small range by varying 

the pressures upstream of the sonic orifices. Figures 11a and b, are tem.nera~ 
ture profiles for hydrazine burning in oxygen-nitrogen atmospheres of different 

' Figure 11a being at a higher value of 
The temperature profile shown in Figure 12 was obtained by reducing the 
oxygen concentration until the only visual evidence of burning was a halo- 
type glow close to the drop surface (Figure 13). Under this condition the 
actual temperatures were lower than those for the other cases but the steep 
gradient at the drop surface is still evident. The change in distance scale 
should be noted. Also it should be noted that the temperature near the drop 
surface with N^H^/air is not much higher than the maximum temperature of 
the halo flame. 

To verify that the temperature profiles obtained for hydrazine were 
functions of the hydrazine itself and not related to the experimental proce- 
dures or equipment, the flame surrounding a heptane drop was also probed. 
Heptane was selected because of its endothermicity . Therefore, if the 
temperature gradient at the surface of the hydrazine drop is due to hydrazine 
decomposition, it would be expected that the profile obtained for the heptane 
drop would not exhibit this steep gradient. Figure 14a is the temperature 
profile obtained for heptane. As can be seen, the temperature rises to a 
peak and then decreases to the droplet boiling temperature. Figure 14b shows 
that two visible "flames" are also present when heptane is burned in nitrogen 
tetroxide. Other investigators , for example, Wharton, Miller, et al, (Ref. 5 
and 6) have also observed a dual reaction zone in premixed n-butane-NO„ 
flames. These results suggest that the two "visible flames" are associated 
with the oxidizer while two real flames are found only for hydrazine. The 
combustion of hydrazine with any oxidizer will be a two flame structure con- 
sisting of decomposition at the droplet surface and the oxidation of the de- 
composition products. 

Concentration Profiles. - The concentration data taken at several 
positions in the flamie are plotted versus reduced radius in Figure 15. It 
is noted that hydrazine was not found in measurable quantities. There is 
some question as to whether this is due to the nonexistence of hydrazine in 
the flame or to hydrazine decomposition in the sampling probe. 

Figure 15 also indicates that the outer flame consists of the decom- 
position of NO^ to NO and O . The mole fractions of NO and O^ decrease 
as the inner flame is approacned. 

The fact that the decrease in NO occurs in the region where the 
tem.perature reaches a maximum is in agreement with the results of 
Reference 5. That investigation demonstrated that the maximum temperature 
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in an n-butane-NO- flame was due to the energy release associated with the 
decomposition of Ml). At the inner flame boundary the oxidizer concentra- 
tions drop sharply to very small amounts. 

Several phenomena which relate to the accuracy of the sampling 
results were observed. It was found that ammonium nitrate was collected 
in the probe throat at the yellow inner flame boundary. This problem 
prevents an accurate description of the inner flame composition being made 
by direct sampling. It Is not apparent whether the ammonium nitrate is 
formed in the probe or Is precipitated from the flame. The formation of the 
ammonium nitrate in the probe would require very fast reaction rates at low 
temperature. Ammonium nitrate is a stable product of the room temperature 
hypergollc reaction between NH_ and NO^/ therefore, it is conceivable that 
it may also be an intermediate or the high temperature reaction. 

The relative amounts of nitrogen were larger than they should have 
been, but this may be due to contamination from the probe purge system. 

CONCLUSIONS 

The facts that, all tests conducted with hydrazine indicate the 
presence of a steep temperature gradient near the drop, and that this property 
is associated only with hydrazine droplets, lead to the conclusion that 
hydrazine decomposition occurs at or very near the drop surface. This 
conclusion appears to be further verified by the concentration profiles obtained 
with the sampling probe. The failure to detect any hydrazine in the flame 
especially in the region near the drop surface, when viewed in the light of 
the temperature profiles, is a further indication of hydrazine decomposition 
at or close to the drop surface. 

The results of these experiments indicate that at low Reynolds number 
the dominant process controlling the hydrazine burning rate is hydrazine 
decomposition, since this mechanism determines the temperature gradient at, 
and therefore the heat transfer to, the droplet surface. The existence of the 
decomposition flame at the drop surface causes hydrazine to have a signifi- 
cantly higher burning rate than thermally stable fuels such as heptane. 
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APPENDIX A 


ANALYTICAL VAPOR PHASE DECOMPOSITION MODEL 
Analytical Model 


Derivation of Equations . - The analysis of the burning of a drop of hy- 
drazine in a nitrogen tetroxide atmosphere has been formulated in a manner 
similar to that described by Coffin and Brokaw t(Ref. 2) ■ in which equilibrium 
chemistry was assumed. As in Reference 2, the problem is set up for a general 
fuel -oxidizer combination with the particular system under examination deter- 
mining the necessary physical properties. 


The equations describing the processes occurring in the region between 
the fuel drop surface and the ambient oxidizer atmosphere are: 

Conservation of Mass 
V 


w. 


w. 


a 


1 ,^ 



l,i 



x= 1 


s 


til 

= Rate of transport of the i chemical species in (moles/sec) 

= Index for atomic species. (e.g.<C=l indicates the 1st 

atomic specie) 

= The number of X atoms in the i^^ specie 


s = Total number of atomic species 

V = Total number of chemical species 


This equation expresses the fact that the transport of the Z atomic species 
through any surface is equal to the sum of the contributions due to the trans- 
port of chemical species containing Z atoms. Since there is no sink in the 
system the net transport of oxidizer atoms is zero. However, there is a net 
production of fuel atoms since the drop is a fuel source. Therefore 


W, 




V\^ = the fuel burning rate in (moles/second) and ^ = the number of Z atoms 
in a fuel molecule. ' 


which describes the conservation of atoms in the systems. 
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Energy Equation 


WA W.H. = -477r® 

f fo 1=1 11 dr 

= The enthalpy of the 1^^ species at temperature T. 

H, = The enthalpy of the fuel at its initial conditions 
fo 

K = Thermal conductivity 
Diffusion Equations 

The equations describing the diffusion of the i species are ob- 
tained from the general diffusion equations (Eqn. 8-1-3 of Ref. 7) by neglecting 
external forces, pressure gradients, and thermal diffusion. In addition it is 
also assumed that each specig^i obeys the ideal gas law and that the binary 
diffusion coefficients of the i specie (D^j) fi\to each of the other species 

are equal, i.e. D = D^. Under these conditions the diffusion equations can be 
put in the form (Ref . 2 ) . * „ i - „ „ 


, 2 P dx. 

dr 


+ — .S, D. + -^ W 

X. J = 1 j dr X, V 


i = 1 , . . .V-l 


R = Gas constant 

. . , .th , , Moles i x 

X = The mole fraction of the i specie ( 77T~ ) 

^ ^°^®fmix) 

For a u component gas there are l^-l independent transport rates. 
Chemistry Equations 


If the transport rates are much slower than the chemical rates the 
burning process can be considered to occur under equilibrium conditions. For 
ay component system there arey-s independent equilibrium relationships 
among the chemical species. ^ „ 

U P. ' 

These relations are of the form ^ 

th '’i 

P. - Partial pressure of the i specie (~^"^x ) 

th ‘ ^ 

- Equilibrium constant for the m reaction. 

th 

„ = Stoichiometric coefficients for the products of the m^" reaction 

i / m 

th 

= stoichiometric coefficients for the reactants of the m reaction, 

X f III 

If the transport and chemical rates are comparable, the relationships 
between the gas components must be determined by finite rate chemistry. 
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Normalization Equation 


V 

S = 1 

i=l 

That is, the sum of the mole fractions is unity. 

In summary, there are s conservation of mass equations, one energy 
equation, v-1 diffusion equations, u -s chemical equations, and a normaliza- 
tion equation for a total of 21 ^ + 1 equations. There are v+1 unknown W's 
(W, i = 1 , . . . . 1 ^ and Wj and u unknown fractions (x, i = 1 , . . . . 1 /) . Thus 

i I 1 

there are 2u+l unknowns and 2v +1 equations and the problem is solvable in 
principle. 

Since there are vW's and u - 1 diffusion equati^c^ns one of the mass con- 
servation equations must be used to determine the v W. The conservation 
equations are 

V 

Ti a. .W. = O'. , W, ^=l,....s 

i,£ 1 f,JL f 

Select i for which ^0 

Designate this value of i as Xl. 

Then 


U-l 

a 

i=l 


. . .W, + a , w = 0 
i,Al 1 V, XI i> 


or 


W = - - 
V Oi 


v-1 

S a. , , W. 

1 , XI 1 


P,Xl i=l 

Substituting for from the diffusion equation and solving for 


p v-1 
477r® P 


W = 
V 


RT 


i=l i,XlV idr^ ^1=1 J d7~^ 


v-1 

D, 


dx. 


i^-l 


X. 


Oi + S Q! 

i^,Xl 1-+1 “i,Xl 


Therefore 


dx 

-W ,=477r^ — : (D. — ^ 

RT L i dr " ’ 


■--I D. , 
x^ j=l j dr 


X. i/-i ax, 

X, v-1 dx, — S a, . (D, - 7 ^+— -tv D, -j^) . 
i. , ^ T- T. X, 1=1 i,Xl 1 dr X J=1 j dr 1 


v-1 


dx, x.v-l dx. 

-i 

j 


X. 


a +7^, a, . , — 
l^,Xl 1=1 i,Xl X 

V 


13 



In order to eliminate all the W's from the energy equation a second mass con- 
servation equation is used to express in terms of the W^'s. Select X for 

which a.^ , ^ 0. Designate this value of ^ X . The l=\ conservation equation 
, f,-t a s 

gives 

W, = — a. . „ W. 

^ “f,X2 ^ ^ ^ 

Substituting into the energy equation gives 

H V 2 

- — a. 0 - S, W.H. =-4 ¥ r K ^ 

i,X2 1 1=1 11 dr 


Rearranging 

S f H. - H, ^ W. ‘V . iA? .., H - H ) W =^irr®K 


The S-2 remaining conservation equations are 

.2-!. OL. = a, . „ w. 

1=1 i,X ^2 i/X 2 1 


^e= 1 ,.. .s 
^ A ^ j '^2 


( 

r 

OL , 

1^1 

^ i,i 

l^-l ^ 

OLj- „ 


_Ul. 
oc. . 

i=l 

f,X2 


Jsl. 

X2 


f,X2 


Equations (1) and (2) are of the form 


l^-l 

S A.W. + BV\^ = C 


Examining the terms in the above equation separately gives for the 
first term 


v-1 


^3 g a.d, + 7 : a.—Ce. d. ^ 

HTTT 1-i 1 11=1 1 1 dr 1=1 ix ^j=i j — -- y 

V •’dr 

S§1 \ X Si °^i,Xl ^i dT" S~x5l "^iirSl “i.XlirSl 

V V V 


X, v-\ 


dx. 

_J_ 

dr 


where u-1 x 

S = i 1 +^i a, . , — i 
V ,X1 1=1 i,Al X 

V 


This relation can be rearranged into the following form 

v-1 

RT 

47Tr ^ P 


■,4i^i^~ ii=i^ i Si \ -r ^§1 h;r O’sSi “i,xiir)}^i 


dx. 
1 

dr 


The second term is 


r,m dx. 

— -- g- BW = ^ a . D — ^ 

4 7Tr^P u Si=l i,Xl i dr 


n*^-l 

^ !£, a. 


x. u-1 dx 
Si=l i,Xlx^j^l j dr 


u-l u-1 x. dx. 

= -r (a. .,+.S, a. ~^D. “T^ 
S 1=1 V i,Xl 1=1 i,Xl x 1 dr 

u 


Therefore Equation (3) becomes 
U-1 _ , U-1 X, u-1 


u-1 


, . _ , , - X. _ O'. . , 0! , . , . dx. 

if, V , 1 -s'' » 1 \ 1^X1 1=1 1/Xl X i \ -n, 1-i 

^ , ; -1 A. +^, A. — - A. — ) , I 1 D. p— I 

1=1 i_ i 1 1=1 id 1=1 IX L TT1 J '* ^ dr _i 

u u V-1 X. 


O' . , O'. . . 1 

t'Al 1/Xl — 

1=1 X 

u 


dT 

For Equation (1) C = - 4Trr^K 


CRT 
47rr^ P 


For Equations (2) C = 0 

Substituting into the above equation gives 


l^-l 


X. 

1 


u-1 u-1 X. u-1 X. ^ O'. O'. , dx 

_i_ ^ 1 , XI 1=1 i,Xl x^ Hi- ^ 


Z-j , i j A.+.^, A. - (B4^ A. ) 
1=1 L 1 1=1 ix^ 1=1 1 x^ ^ 


l^-l X, 

1=1 1,X1 x^ 


J J i dT 


= C 


(4) 
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where 


D. 

6 = — 
i KRT 


For Equation (1) a, ^ _ 

A, 

' ^ f,X2 

%,X2 


fo - 


B _ H, - H. 

f,X2 fo 1 


■N 

J 


C' = p 

For Equations (2) 


A. 


B 


= (a, 


i,Si a. 


= (a, 



“i,X2^ 

i = 1 
\ 

°^f,X2 


\,X2^ 


“f.X2 



C =0 

The chemical equilibrium relations are 


€. 


K 


? Pi 


i,m 


m 


n p. i,m 
1 1 


Taking the temperature derivative of this equation gives 


« o ciK 

e V B. V p. E 

^ n p = K ^ n p. + ^p, dT 

dT i ^i m dT X 1 


V 


T^i 


u 


dZ, 'i,m Vt 

? Pi ' ? Pi 'l,m Pi 


-1 ^ 
dT 


Substituting and rearranging yields 


1/ ... (e. -B. ) 

^ : e II pK. - K B. 

i=ltLi/mj j mi, m_i x 


n “^i 1 


dK 


m 


dT J dT 
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(5) 

m = 1 , . . .1^ -S 


Where the relationship between the partial pressure and the mole fraction, 

X, = P,/Pr has been used. 

1 1 


It should be noted that Equations (5) provide i/-S relations between the 
X 's. While these equations were derived for equilibrium chemistry and other 
chemical assumptions, e. g. finite rate kinetics, or arbitrary specification of 
X. vs T, can be made, providing a set oft'-S independent relationships are 


obtained . 


The normalization equation can also be differentiated to give 


V 

E 

i=l 



( 6 ) 


Equations (4), (5), and (6) are respectively S-1, I'-S, 1, equations for a 
total of V equations for the v unknown x^'s. 

The above equations can be numerically Integrated to give x^ vs T for 
all the species. The mass conservation equations can then be used to obtain 
the spatial specie distribution from the drop to the ambient atmosphere. 


Since the experimental programs described above indicated that the 
hydrazine burning rate was controlled by decomposition close to the drop surface 
the usefulness and applicability of an analytical model of this type toward the 
prediction of the burning rate was uncertain and a computer code was not 
developed. 
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APPENDIX B 


PHYSICAL PROPERTIES 

For use in the analytical hydrazine decomposition model, enthalpies, 
diffusion coefficients, and thermal conductivities are needed for the following 
chemical species: N 2 ^ ^ 2 ©, NO, NO 2 / O 2 / O, N, NH^NO^/ 

and N 2 H^. 

Enthalpies . - The enthalpies of H 2 / NH^/ N 2 , H 2 O, NO, NO 2 / O 2 / 


O, and N are given in Reference ( 8 ) 


from 0°K to 6000° K. The enthalpies 


of N„H, and NH.NO_ are not given but can be determined from the same 
2 4 4 3 

NASA-Lewis program. However, apparently the program shows that N 2 H 4 does 
not exist as such at any temperature since equilibrium calculations are involved, 

Diffusion Coefficients. - Two suitable expressions are available for 
determining the coefficient of diffusion in a binary mixture. 

Hirschfelder, Curtiss, and Bird, (Ref. 4 ) give the analytical expression 

^T^CM +M )/2M M 

D,„ = 0.0026280 

pCTiz « 12 ^^12 ' 


where 


= diffusion coefficient 

= pressure 

C7 (T 

1+2 


cr^,a 2 = molecular diameters of species 1 and 2 


= temperature 


Mf = molecular weights of species 1 and 2 

(1 1 ) * 

f 2 j^ 2 ^ “ integral parameter, approximately equal to 1 . 

★ 

Tj ^2 “ overall reduced temperature between the two species 

Another suitable expression for estimating the diffusion coefficient between 
species 1 and 2 is the empirical Gilliland (Ref. 9 ) expression 
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D = 0.0043 


Ml 

where 

D = diffusivity, cm*/sec 

T = temperature, °K 

p = pressure atm 

M 1 M 2 = molecular weights of gases 

V^,Vg = molecular volumes at the normal boiling points, cc/g-mole 

The molecular volume can be calculated from the molecular weight and 
the density of the substance at the boiling point. If the density at the boiling 
point is unknown, the molecular volume can be calculated from the sum of the 
atomic volumes. The atomic volumes are given in Reference 10. 

Thus, if the molecular diameters are given (some are given in Ref. 7) 
the more exact Hirschfelder, Curtiss, and Bird equation can be used. If the 
molecular diameters are unknown, the Gilliland equation can be used. Either 
equation will give a reasonable estimate of the diffusion coefficient in a 
binary mixture. 

Since a binary gas mixture is assumed in use of the equations, in the 
actual multicomponent case the diffusion coefficient can be obtained from a 
mass^jwelghted average of the properties of the remaining species into which 
the i specie is diffusing. If the specie diffuses into a mixture containing a 
large amount of one component with respect to the others in the mixture (l.e. 
approaching an overall binary mixture between the specie and the remaining 
material) the calculation will be reasonably accurate. The error involved if the 
overall mixture (specie diffusing plus the remaining materials) is not a binary 
mixture is not known at this time. It should be mentioned, that in a binary 
gas mixture, the dependence of the diffusion coefficient on composition is 
only slight. 

Thermal Conductivities . - Since only limited experimental thermal 
conductivity data is available, it is necessary to use a theoretical expression 
to determine the thermal conductivity. 

The thermal conductivity of a substance can be reasonably estimated 
from an analytical expression given by Hirschfelder, Curtiss, and Bird, 

(Ref. 7). 



19 


k = 1989.1 X 10 


j^(2.2)*(t*) ^15 


^ 1 

R 3 y 


k = thermal conductivity 

T = temperature 

M = molecular weight 

(T = molecular diameter 

(2 2 )* 

Q = Integral parameter# approximately equal to 1 . 

C = heat capacity 

P 

R = universal gas constant 

T* = reduced temperature 

As used in the model, the thermal conductivity, k, is Included because 
of a comparison with heat conduction. Therefore, k represents an overall 
constant value. In this analysis, the thermal conductivity will be determined 
from a mass weighted average of the thermal conductivities of the respective 
species. The closer one specie comes to being present in a large amount, the 
more accurate the use of the mass weighted average thermal conductivity will 
be. 

Whenever possible experimentally determined values of the thermal con- 
ductivity of the substance should be used. 

High temperature thermal conductivity data (3000-4700°F) is available 
in Reference 11. Figure 16 shows experimental results obtained from Reference 
11 in addition to results from Reference 12 (indicated by the broken line). 


The following low temperature experimental thermal conductivity (k) data 
is available for all of the species except O, N, NH^NO^/ and N2H^* 


Temperature , 

calories 
^ cm sec°C 

Refer 

- 57.6°C 

3.82 

13 

0 

5.135 

13 

100 

7.09 

13 

-191.4 

1.829 

13 

- 78.4 

4.305 

13 

0 

5.68 

13 

100 

7.18 

13 
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Temperature 

calories 
cm sec C 

Reference 

^2 

-252. 2°C 

3.22 

13 


- 78.4 

30.65 

13 


0 

39.60 

13 


100 

49.94 

13 

H„0(g) 

46 

4.580 

13 


100 

5.510 

13 

NO (g) 

- 71.4 

4.160 

13 


0 

5.55 

13 

N02(g) 

55 

8.88 

13 

°2 

-191.4 

1.721 

13 


- 78.4 

4.292 

13 


0 

5.70 

13 


100 

7.427 

13 


80°K 

1.701 

14 


90 

1.930 

14 


100 

2.159 

14 


110 

2.387 

14 


120 

2.614 

14 


130 

2.840 

14 


140 

3.064 

14 


150 

3.287 

14 


160 

3.508 

14 


170 

3.728 

14 


180 

3.946 

14 


190 

4.162 

14 


200 

4.375 

14 


210 

4.584 

14 


220 

4.790 

14 


230 

4.993 

14 


240 

5.194 

14 


250 

5.392 

14 


260 

5.586 

14 


270 

5.780 

14 


273.1 

5.839 

14 


280 

5.970 

14 


290 

6.159 

14 


293.1 

6.218 

14 


298. 1 

6.314 

14 


300 

6.350 

14 
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Temperature 

, calories 
cm sec C 

Refer 

310°K 

6.547 

14 

320 

6.748 

14 

330 

6.954 

14 

340 

7.164 

14 

350 

7.378 

14 

360 

7.594 

14 

370 

7.812 

14 

380 

8.033 

14 


Thermal Decomposition of Hydrazine. - In connection with the model/ 
information on the thermal decomposition of hydrazine is desired. Although 
there are discrepancies in the literature/ the following information is considered 
significant. 

Thomas (Ref. 15 ) shows that hydrazine decomposes according to two 
consecutive overall reactions at 20 atm total pressures 

3H2H^ -♦ 4NH2 + N2 + SO.lSKcal ( I) 

4NH2 2N2 + 6H2 - 44.00Kcal (II) 

The overall decomposition reaction is 


3N2H^ -* 4(1 -x)NH 2 + (l+2x)N2 + 6XH2 + (80.15-44.00x)Kcal 
where x = the fraction of NH^ decomposed. 

Ammonia (NH_) is only negligibly dissociated for temperatures below 400°K. 

As the temperature is raised above 400°K, NH^ decomposes with the decomposition 
becoming complete at 800°K. Experimental results show that the decomposition 
reaction of NH (II) is slower than reaction I. Hence/ if the decomposition (I) occurs 
in the order of milliseconds , it then might be reasonable to assume that only 
part of the NH^ could decompose. 

Lucien (Ref. 16 ) carried out thermal decomposition studies of hydrazine 
in the ranges of 175° to 250° C and 300 to 430 psl. He found that the decom- 
position went according to 

3N2H^ -♦ 4(1-x)NH 2 + (l+ 2 x)N 2 + 6 XH 2 
where x = 0.017 at 222° C and 300 psi 
= 0.034 at 250° C and 430 psl 
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The calculated activation energy was about 73Kcal per mole which is greater 
by a factor of two than other reported values. The rate of decomposition 
decreased with pressure over the range studied. 

Eberstein and Glassman (Ref.l7) studied hydrazine decomposition. The 

rate constant for hydrazine was found to be 

1 n 33 -1 

k = 10 exp (-36,170/RT)sec 

The stoichiometry observed was 

N2H^ 0.9 NH^ + O. 5 N 2 + O. 6 H 2 

The disagreements in the literature may indicate that the hydrazine 
decomposition mechanism changes with pressure and temperature. 
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■ FIGURE 1. HYDRAZINE / NITROGEN TETROXIDE 


26 







Ih 

0) 

N 

i-, 

O 


10 

> 


O 


CP 

c 

•<— < 

a 

S 

ro 

CO 


<D 

> 

*rH 

i-4 

Q 

0) 

s: 

o 

S-I 

Oh 


CD 

CP 

C 

CO 


CD 

> 

Q 



27 


FIGURE 2. HYDRAZINE DROPLET BURNING SETUP. 
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FIGURE 9. HYDRAZINE/NITROGEN TETROXIDE FLAME 
TEMPERATURE PROFILE 
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FIGURE 12. HYDRAZINE/LOW OXIDIZER CONCENTRATION 
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FIGURE 13. HYDRAZINE/LOW OXIDIZER CONCENTRATION 
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